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The increasing sensitivity of cosmic microwave background (CMB) missions will require to sig- 
nificantly improve the accuracy in the subtraction of the various sources of Galactic foreground, 
from the most relevant components (synchrotron, dust and free-free emission) to those usually 
considered of minor relevance in CMB experiments. With respect to other Galactic diffuse com- 
ponents, the Zodiacal Light Emission (ZLE) is peculiar, depending not only on the observing 
direction but also on the location of the observer within the Solar System: ZLE behaves then 
as a large scale, time-dependent foreground. Starting from the existing far-infrared ZLE mod- 
els, we discuss the impact of ZLE contribution in CMB maps and the level of contamination in 
time ordered data and maps expected from the forthcoming PLANCK space mission as well as the 
Planck capability to increase our knowledge of the ZLE properties. 
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1. Introduction 

The increasing sensitivity of cosmic microwave background (CMB) missions will require a 
careful subtraction of any source of weak "Galactic" foreground, asking to consider in CMB ex- 
periments other sources of Galactic foregrounds over the most relevant ones, i.e. synchrotron, dust, 
and free-free emission. In this context, the Zodiacal Light Emission (ZLE) due to the thermal emis- 
sion from the cloud of Interplanetary Dust Particles (IDPs) permeating the Solar System represents 
a significant component. At frequencies rs 10 THz the ZLE dominates the sky emissivity, and in 
preparing templates of Galactic emission from observations at these frequencies the contamination 
of ZLE has to be accurately accounted for, since the ability to model and remove the ZLE will 
largely fix the final accuracy of the Galactic templates. At frequencies below rs 1 THz the ZLE is 
subdominant compared to the Galaxy, but its surface brightness is still significant particularly in 
regions where the Galaxy emission is weak [fil 01. It has also to be considered that the ZLE does 
not depend only on the instrument pointing direction, P, but also on the position of the observer, 
Rp, within the Solar System. Then, the ZLE behaves as a time-dependent foreground and, when 
not properly removed, introduces subtle systematic effects. 

2. Model and extrapolations 

As for the Galactic dust contamination, the modelling of the ZLE has to be based on far- 
infrared observations, at least in what regard the geometrical aspects of the IDPs distribution. Key 
data for the ZLE below 300 jiim have been obtained by IRAS [g], COBE [|, g] and ISO |@, |]. The 
starting point for this analysis is the COBE/DIRBE model for the ZLE which describes the expected 
local 3D emissivity within the IDPs complex [Q]. In particular, among the various components in 
which the IDPs are distributed we refer here to the dominant Smooth component which accounts 
for more than 90% of the ZLE. Then, for frequencies f ^ I THz the total brightness of the ZLE 
integrated along a given line-of-sight is 

/ / (P,Rp)=£ / Z / (P,R P ), (2.1) 

where Z/(P,Rp) gives the spatial dependence and the Emissivity Factor Ef is a correction with 
respect to a pure blackbody emission law, related to the composition and size distribution of dust 
grains. Following [Q] it is assumed £y=lfor/=12 THz. The spatial dependence is given by the 
integral along the line-of-sight 

Z / (P,Rp)=/ dsN(R P + sF)B f {T(R P + sF)), (2.2) 
Jo 

where s is the distance from the observer along P, Bf(T) the blackbody brightness, N(r) the dust 
density at a given location in the Solar System (assuming A^(r) = for |r| > 5.2 AU), T(r) the 
local dust temperature, assumed to scale as |rj~ - 467 . 

In long duration CMB experiments the sky is observed while the observer is orbiting around 
the Sun and consequently the ZLE will show seasonal modulations at the level of 5 — 15% [Q]. 
Changes in the scanning strategy or observation epoch will result in different realizations of ZLE 
sky maps even when the same set of pointing directions is taken. To condensate into a statical map 
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Table 1: Predicted values of Ef and of the ZLE at the highest PLANCK frequencies. The table reports 
minimum, mean and maximum ZLE within a w 50% uncertainty for a one year mission. 



this dynamical information we exploit the cylindrical symmetry of COBE distribution of the IDPs 
Smooth component. We develop a serie expansion of Z/-(P,Rp) about an averaged orbit in the 
IDPs cloud reference frame [ph. This is equivalent to calculate a kernel map for Zy(P,Rp) which is 
good for a given "nominal" mission together with coefficients to be used in computing variations 
of the map for a range of possible variations in the mission. This method allows at the same time 
the generation of data streams of ZLE signals for a given list of spacecraft positions and pointing 
directions, and the generation of specialized time averaged maps (e.g. yearly averaged maps) for a 
specific mission, orbit and scanning strategy 1 . 

The extrapolation of Ef at frequencies below 1 THz is a more delicate problem. COBE/DIRBE 
measures extend down to / = 1.2 THz, which fixes the lowest frequency for which the COBE 
model provides values of Ef [Qj. However COBE/FIRAS provides measures of ZLE averaged 
over the sky and over one year down-to / « 3 x 10 2 GHz [jl|] but with a not high S/N ratio. We 
then compare the COBE/DIRBE measures with simulated yearly and full-sky averaged values of 
Z/(P,Rp) in order to evaluate Ef. In this way we can obtain numerical estimates for these parame- 
ters [§]. Similar results are also obtained directly extrapolating Ef values derived by COBE/DIRBE 
at / < 12 THz down to the required frequencies. Tab. [I] summarizes the results of these calculations 
giving the minimal and maximal yearly averaged levels of contamination at frequencies relevant 
for the Planck mission. 

The left frame of Fig. [I] represents a typical realization of a ZLE signal at 857 GHz (together 
with secondary components and uncertainties) compared to the Galaxy. The right frame shows the 
relative contribution of ZLE brighness and instrumental noise respect to the Galaxy. In the plot 
the ZLE and the noise contributions are averaged over a circular band of about 85° of radius and 
drawn around an axis of given ecliptical longitude. The same is done for the Galaxy. The red 
band is the noise, the white line the expected ZLE over Galaxy averaged ratio, the blue band the 
±lrj band, the yellow line the maximum ratio. The plotted values are ordered as function of the 
ecliptical longitude of spin axis vector. The data are calculated for patches of 1° in radius. The 
noise is referred to a 14 month mission (2 sky surveys). 



Tables of serie expansions of ZLE for all of the PLANCK frequencies, the related software in IDL, and documen- 
tation may be required to the authors. 
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Figure 1: Left panel: simulated data stream of surface brightnesses (MJy/sr) expected at 857 GHz for the 
ZLE - smooth component (red), secondary components (cyan) and uncertainties in the prediction (bands) - 
the Galaxy (green), and the sum of ZLE and Galaxy (blue). Also signals from ZLE (red-dashed) as after 
summation of the Galaxy (green dashed) are shown. The ordinate is the phase along the scan circle. Right 
panel: relative averaged contribution of ZLE, and noise compared to the averaged Galaxy and computed 
over circular bands, as a function of the ecliptical longitude of the spin axis vector. 

3. Detection and removal of the ZLE 

It is evident that current information does not allow to reach an accuracy better then as 20% 
in the ZLE removal; this is due in particular to the uncertainties on Ef below 1 THz. New, direct 
measures are necessary to improve the situation. The ability of PLANCK to detect the ZLE signal 
improving the accuracy in Ef determination has been analyzed in [§] assuming that the COBE 
model properly represents the spatial distribution of the ZLE leaving as a free parameter Ef. Two 
methods have been considered, the first based on the comparison of a template map for the Galactic 
emission with a spatial template for the ZLE calculated for given mission orbit and scanning strat- 
egy, the second method is based on the comparison of observations of the same regions of sky taken 
at different epochs then exploiting the seasonal dependence in the ZLE observation. In both cases, 
since the ZLE varies over scales of « 10°, one can consider template maps and observations at 
resolutions of 1° — 2° (this alleviate also the contribution of local features, as weak point sources 
or regions with peculiar frequency dependencies). In this kinds of analysis, it is important to apply 
cuts to the data excluding regions where the Galaxy is very bright. Simulations shows a typical 
absolute RMS uncertainty on Ef determination induced by the limited instrumental sensitivity of 
~ 1(T 3 , 2.1 x ~ 1(T 3 and 2.6 x ~ 1(T 3 at 857 GHz, 545 GHz and 353 GHz, respectively. For 
typical expected values of E f (« 0.65, 0.26, 0.11 for 857 GHz, 545 GHz, 353 GHz) the Planck 
sensitivity will allow an Ef recovery at 0. 15%, 0.8% and 2.4% (la) accuracy at 857 GHz, 545 GHz 
and 353 GHz, respectively. Of the most relevant systematic effects, pointing and sampling uncer- 
tainty, aberration of light, Doppler shift and relative calibration uncertainty, only the last one is 
found to be really critical, while the errors introduced by the other ones are found to be signifi- 
cantly below the noise. Then, relative calibration could ultimately determine the final accuracy in 
the ZLE extraction from PLANCK data. For a relative calibration RMS error of ~ 1% (0.1%) on 
patches of 2° radius, we find an absolute RMS error on E f of ~ 0.01 - 0.04 (~ 0.001 - 0.004) with 
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only a weak dependence on the frequency in the range ~ 300 — 900 GHz, corresponding to relative 
errors on E f ~ 4%, 10%, 23% (~ 0.4%, 1%, 2%), respectively at 857 GHz, 545 GHz, 353 GHz 
for the most likely Ef values expected on the basis of COBE/FIRAS data. 
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